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HIGHLIGHTS 


• We review experimental data from a VI60 engine developed for cogeneration. 
. We also investigate the V161 solar engine. 

• The possible margin of improvement is evaluated for each subsystem. 

• The procedure is based on similarity models and thermodynamic models. 

. The procedure may be of general interest for other prototypes. 
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The development of systems based on Stirling machines is limited by the lack of data about the perfor¬ 
mance of the various subsystems that are located between the input and output power sections. The 
measurement of some of the variables used to characterise these internal subsystems presents difficul¬ 
ties, particularly in the working gas circuit and the drive mechanism, which causes experimental reports 
to rarely be comprehensive enough for analysing the whole performance of the machine. In this article, 
we review experimental data from a VI60 engine developed for cogeneration to evaluate the general 
validity; we also investigate one of the most successful prototypes used in dish-Stirling systems, the 
VI61 engine, for which a seemingly small mechanical efficiency value has been recently predicted. The 
procedure described in this article allows the possible margin of improvement to be evaluated for each 
subsystem. The procedure is based on similarity models, which have been previously developed through 
experimental data from very different prototypes. Thermodynamic models for the gas circuit are also 
considered. Deduced characteristic curves show that both prototypes have an advanced degree of devel¬ 
opment as evidenced by relatively high efficiencies for each subsystem. The analyses are examples that 
demonstrate the qualities of dimensionless numbers in representing physical phenomena with maxi¬ 
mum generality and physical meaning. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

Features of the Stirling engine that are of greater interest for 
current applications are most likely the ability to adapt to various 
energy sources and to provide not only mechanical or electrical en¬ 
ergy but also deliver thermal energy to its cold sink with low levels 
of noise and chemical pollution. 

Systems based on Stirling engines contain various subsystems 
connected between the input and output power sections. Power 
is supplied by a heat exchanger adapted to the characteristics of 
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the energy source. The useful power is frequently obtained by an 
electrical generator driven by the engine mechanism and, in cogen¬ 
eration applications, through a heat exchanger powered by the en¬ 
gine’s cold source. 

To characterise each subsystem, a number of variables must be 
measured, preferably by direct measurement under real operating 
conditions. However, the measurement of some variables presents 
difficulties, particularly in the working gas circuit and the drive 
mechanism because these subsystems are not directly connected 
to the input and output power sections. Consequently, it is rare 
that experimental reports are comprehensive enough, which hin¬ 
ders the validation of simulation procedures and the development 
of applications. 
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Notation 



Acx 

cross-sectional area of space x (m 2 ) 

v sw 

engine swept volume = V max - V min (m 3 ) 

^WX 

wetted surface of space x (m 2 ) 

a R 

regenerator thermal diffusivity (m 2 /s) 

Cr 

specific heat of regenerator material (J/kg I<) 

&CX 

cross-sectional area ratio of space x = Acx/V 2 ^ 

G w 

water flow at secondary cooling circuit (m 3 /s) 

Palt 

alternator efficiency = P e /P B 

h 

heat loss factor of control volume 

Pb 

brake effkiency= Pb/Qe 

Lc 

cooler tube length (m) 

Pcomb 

combustion subsystem efficiency= Q.e/P/ n 

Le 

heater tube length (m) 

>1e 

total electrical efficiency = P e /P, n 

Lr 

regenerator length (m) 

le 

electrical efficiency based on the heat power input to 

Nma 

characteristic Mach number = n s V]^ / ^RT wC 


heater = P e /Q_E 

N P 

characteristic pressure number = N S cN M a = Pm^lw / 

Pind 

indicated efficiency =P lnc //Q.£ 


(HVRT„c) 

Imec 

mechanical efficiency = P B /P !n d 

Nsc 

characteristic Stirling number = p m /(pn s ) 

Ur 

thermal efficiency at primary cooling circuit =P f /P,„ 

N T cr 

characteristic regenerator thermal capacity number 

nut 

thermal efficiency at secondary cooling circuit = P ut /P jn 


= PR c RTwc/P m 

<s> 

dimensionless coefficient of linear losses of indicated 

K 

characteristic regenerator thermal diffusivity number 


power 


= °W(V’4 3 '/KT w d 

y 

adiabatic coefficient 

n s 

rotational frequency (rev/s) 

K 

swept volume ratio = V c /\4 

Pb 

brake power (W) 

hh 

dimensionless geometrical parameters, including those 

Pe 

electrical power output (W) 


characteristic of the drive mechanism 

Pin 

fuel equivalent thermal power input (W) 

Khx 

hydraulic radius ratio of space x = r^/V]^ 

Pi nd 

indicated power (W) 

P-dx 

dimensionless dead volume of space x = Vdx/Vsw 

P 

1 mec 

mechanical power losses (W) 

Kv 

regenerator volumetric porosity 

P t 

thermal power output at primary cooling circuit (W) 

<F 

dimensionless coefficient of quadratic losses of indi- 

Pur 

useful thermal power output at secondary cooling 


cated power 


circuit (W) 

P 

working fluid density (kg/m 3 ) 

Pm 

working gas mean pressure (Pa) 

PR 

regenerator material density (kg/m 3 ) 

0.E 

heat power input to heater (W) 

T 

Temperature ratio = T wC /T w£ 

R 

specific gas constant (J/kg K) 

Cind 

dimensionless indicated power = P !n d/(PmV sw n s ) 

Oix 

hydraulic radius of space x (m) 

s mec 

dimensionless mechanical power losses = PmeJiPmVsw- 

T a 

environmental air temperature (K) 


n s ) 

Tc 

cooling fluid temperature at cooler outflow (K) 

Co 

dimensionless quasi-static indicated work per cycle 

Tew 

cold water inlet temperature at secondary cooling 
circuit (K) 

Subscripts 

Te 

maximum working gas temperature at heater head (K) 

C 

cooler 

Thw 

hot water outlet temperature at secondary cooling 

cc 

compression cylinder 


circuit (K) 

ec 

expansion cylinder 

Twx 

wall temperature at space x (K) 

E 

heater 

v dx 

dead volume of space x (m 3 ) 

R 

regenerator 

V c 

compression cylinder swept volume (m 3 ) 

max 

corresponding to the operating point of maximum 

V E 

expansion cylinder swept volume (m 3 ) 


indicated power 


The University of Oviedo and the IK4-TEKNIKER research centre 
have collaborated since 2007 on the development of Stirling en¬ 
gines capable of being used for solar energy conversion [ 1 ] and mi¬ 
cro-combined heat and power (CHP) [2], Since the beginning of 
this collaboration, it was considered interesting to analyse the per¬ 
formance of engines potentially useful for these applications using 
dynamic similarity models previously developed from experimen¬ 
tal data on indicated power and mechanical losses observed in 
benchmark prototypes of very diverse characteristics: General Mo¬ 
tors GPU-3 engine operating with helium and hydrogen, Philips 
M102C engine operating with air, United Stirling P40 engine oper¬ 
ating with hydrogen, helium and nitrogen, Kockums V4275R Mark 
II engine operating with helium, Thermomotor 360-15 engine 
operating with hydrogen, Yamanokami 1 and Yamanokami 2 
low-temperature differential (LTD) engines operating with air, Eco- 
boy SCM81 engine operating with helium and nitrogen, Mitsubishi 
and Daihatsu large-bored engine operating with helium, Volund 
SM-1 engine operating with helium, Kolin’s and Gros’ LTD engines, 
Stirling’s engine of 1818, etc. [3-5], Such models have proven to be 
useful for design using scaling techniques [6,7], but have also 
provided analysis criteria, for example, to delimit the power pro¬ 
duced by the original Stirling air engine of 1818 [8], to evaluate 


predictions of simulation programs [9-11], to carry out compara¬ 
tive studies among empirical models [12,13], and to define operat¬ 
ing ranges in control systems [14], 

From a practical standpoint, experimental tests on Stirling en¬ 
gines are considered sufficient when at least the electrical and 
thermal powers are measured in the output section. Generally, 
such tests are performed for various mean pressure values of the 
working gas, which is often the main control variable. However, 
to evaluate the improvement margins of the subsystems it is desir¬ 
able to obtain their individual efficiencies as a function of as many 
parameters and operating variables as possible. The rotational fre¬ 
quency of the engine is an essential performance variable, but it is 
generally only varied across a small range in grid-connected sys¬ 
tems, which hinders analyses. In this article, the internal subsys¬ 
tems of a Stirling unit are characterised from measurements 
made at the input and output sections. The experimental data 
presented refer mainly to a VI60 Stirling engine developed for 
cogeneration, but an important objective of the article is to draw 
conclusions applicable to one of the most successful prototypes 
used in dish-Stirling systems, the V161 engine, which had its 
mechanical efficiency questioned by a recent study [15], To the 
best of our knowledge there is not any publication about other 
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systems for which the five major subsystems (combustion subsys¬ 
tem, gas circuit, drive mechanism, electrical subsystem and cooling 
subsystem) have been characterized experimentally. Using dimen¬ 
sionless variables facilitates the procedure, which may also be of 
general interest for other prototypes. 

2. Experimental data 

The VI60 Stirling engine has evolved since its original develop¬ 
ment by United Stirling AB, Sweden, in the late 1960s. Develop¬ 
ment was later transferred to Stirling Power System, USA, which 
manufactured a number of versions before reaching commercial 
quality with the “G” series in 1989. The version known as V160- 
DMA is a Stirling cogeneration system that belongs to the “E” series 
[16], Schlaich, Bergmann & Partner, Germany, bought a license 
from Stirling Power Systems and the manufacturing license was 
then transferred to SOLO Kleinmotoren GmbH, which produced 
commercial versions of the VI60 engine in recent decades. The 
VI61 Stirling engine is a VI60 unit with modifications intended 
to reduce production costs, which is adaptable for solar applica¬ 
tions based on the dish-Stirling concept [17], Current models are 
marketed by Clean Energy, Sweden. 

Experimental tests of the V160-DMA unit operating with natu¬ 
ral gas as fuel and helium as the working fluid have been published 
[18], In these tests, the engine subsystems (Fig. 1) were monitored 
to obtain both electrical and thermal output power, as well as the 
corresponding efficiencies. 

The measured data were tabulated in series corresponding to 6 
values of the water flow at the secondary cooling circuit, decreas¬ 
ing from G w = 1.95 m 3 /h (series “a”) to G w = 0.45 m 3 /h (series “f”). 
Each series contains data obtained for 11 values of the main control 
parameter, namely the mean working gas pressure, which was var¬ 
ied from p m = 3.7 MPa to p m = 12.5 MPa. An additional engine con¬ 
trol is obtained through the regulation of the cooling outflow 
temperature T c , which is somewhat lower than the cooler wall 
temperature, while the maximum temperature of the working 
gas is maintained at a near-constant of 625 ± 1 °C. The operation 
of the asynchronous electrical machine causes the rotational fre¬ 
quency to stay near-constant as well, between the idle value of 


1503 rpm at p m = 3.7 MPa to 1518 rpm at p m = 12.5 MPa. Fig. 2 
shows that, for each series, electrical and thermal power measure¬ 
ments are linearly proportional to p m at both the input or output 
sections. 

In addition to power and efficiency data, the tables also include 
measurements of the temperatures T c , and T hw at the primary 
and secondary cooling circuits for each series. Fig. 3 shows that 
these temperatures have a high degree of linearity as a function 
of p m . From Figs. 2 and 3, it is deduced that varying the water flow 
G w at the secondary cooling circuit does not provide a wide range 
of power control, because the trend lines lead to similar power val¬ 
ues for each mean pressure value. This observation is in agreement 
with results obtained by Gheith et al. [19], who observed that flow 
rate has less significant effect in power brake performance than 
heating temperature and charge pressure. 

While a detailed series of environmental air temperature mea¬ 
surements was not provided, the authors reported that they oscil¬ 
lated between 19 and 26 °C. 

3. Analysis of the external combustion subsystem 

The performance analysis that will be described in the fol¬ 
lowing sections requires experimental knowledge of the heat 
power absorbed by the engine heater. In steady operation, the 
combustion subsystem efficiency i] CO mb can be obtained by 
means of the balance of power transferred through the control 
volume (Fig. 1) that includes the subsystems formed by the Stir¬ 
ling engine, the primary cooling circuit and the alternator, that 
is to say: 

>lcomb = 'le + rit + k L ( 1 ) 

where the k L factor has been introduced to consider heat losses 
through the control volume’s casing. 

This k L factor is not available from the experimental data, but it 
can be estimated. Table 1 shows the mean combustion efficiency 
values, their standard deviations, and the mean k L factors, which 
have been obtained for each series, assuming that losses to envi¬ 
ronmental air can be represented by a heat transfer coefficient of 
10W/(m 2 °C), that the control volume casing has a surface area 


p -0 

COMBUSTION 

SUBSYSTEM 


(1 -l\coJP,n 

3 


_Q QE=n mmb Pi„ 





Fig. 1. Power flow between the subsystems. 
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Fig. 2. Variation of input and output powers as a function of the mean pressure. 
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Fig. 3. Variation of temperatures in the primary and secondary cooling circuits as a function of the mean pressure. 


of 0.75 m 2 [20] at a temperature 5 °C higher than the cooling fluid 
temperature T c at the cooler outflow, and that the environmental 
air temperature is 20 °C. These assumptions lead to mean r) comb val¬ 
ues of approximately 90%, with relatively low standard deviations 


for each series. It is deduced that iy com i, is almost independent of the 
mean pressure and therefore the input and output powers; thus, 
the heat power Q £ absorbed by the engine heater is also linearly 
proportional to p m for each series (Fig. 4). 
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Table 1 

Mean values of the combustion subsystem efficiency and k L factor. 



Series “a” 

Series “b” 

Series “c” 

Series “d” 

Series “e” 

Series “f ’ 

Mean efficiency 

Mean k L factor 

0.89 ± 0.01 

0.014 ±0.001 

0.92 ± 0.02 

0.013 ±0.001 

0.92 ± 0.01 

0.013 ±0.001 

0.90 ± 0.02 

0.015 ±0.001 

0.90 ± 0.01 

0.015 ±0.001 

0.89 ± 0.03 

0.017 ±0.001 



Fig. 4. Variation of heat power absorbed by the engine heater as a function of the mean pressure. 


4. Analysis of the internal subsystems 

For each series, it can be written that: 

We Walt Wmec Wind Wcomb 

Thus, once k L has been estimated, the experimental data will al¬ 
low the electrical efficiency to be deduced based on the heat power 
absorbed by the engine heater, i.e., the product >] ind ij mec i] all for each 
value of the mean cycle pressure p m . These data can be used as a 
starting point for the analysis of the subsystems’ performance, 
with the goal of obtaining the subsystem efficiencies as a function 
of the main engine operating variables, i.e., mean pressure, rota¬ 
tional engine frequency and temperatures. 

Series “e” has been selected as the case study. This series’ data 
are listed in Table 2, with an additional four columns for the origi¬ 
nal measurements reported by Lista [18], The values of n s are esti¬ 
mates deduced from the aforementioned operational range of the 
rotational frequency, while the last three columns are a conse¬ 
quence of Eq. (1). 


Fig. 5 shows evidence that the VI60 system reaches tj* values 
which are slightly higher than 26% at a mean pressure ranging 
from 10 to 12.5 MPa. 

4.1. Semi-empirical models for the analysis of subsystem performance 

Organ [21 ] and Prieto et al. [3,22] independently introduced the 
following functional relationship to express the dimensionless 
indicated power of kinematic Stirling engines: 

Cind fir.K.A]..... j-k ■ Pdi't: ■ I L dF ■ hdR ■ hdC - J cE, -'Ac ■ 'liF. ■ 'hR ■ '-hC ■ Ntcr,JV p ,Nma) 

(3) 

Prieto et al. [4] proposed the following approach that explicitly 
states the influence of the engine speed on the dimensionless indi¬ 
cated power: 

C = Co - W>N ma - WN 2 ma (4) 

In this equation, Co denotes the quasi-static dimensionless indicated 
work per cycle that represents the theoretical prediction, assuming 


Table 2 

Series “e” experimental data (G w = 0.75 m 3 /h, T h = 625 °C). 


P m (MPa) 

n s (rpm) 

T c (°C) 

Thw (°C) 

Tew (°c) 

Put (W) 

Pel W) 

Pin (W) 

r]t (-) 

Pt (W) 

Qe(W) 

klPin (W) 

t lcomb ( ) 

3.7 

1505.0 

34.6 

32.2 

26.0 

5310 

1250 

8660 

0.72 

6235 

7629 

144 

0.88 

4.0 

1505.4 

35.8 

33.3 

26.7 

5720 

1310 

9040 

0.74 

6690 

8153 

153 

0.90 

5.0 

1506.9 

40.0 

37.2 

29.3 

6900 

1930 

11,180 

0.70 

7826 

9941 

185 

0.89 

6.0 

1508.4 

43.4 

40.3 

31.2 

8050 

2500 

13,140 

0.69 

9067 

11,777 

210 

0.90 

7.0 

1509.9 

47.1 

43.6 

32.7 

9380 

3160 

15,200 

0.68 

10,336 

13,734 

238 

0.90 

8.0 

1511.4 

50.1 

46.3 

34.3 

10,500 

3750 

17,190 

0.67 

11,517 

15,528 

260 

0.90 

9.0 

1512.8 

53.5 

49.5 

36.4 

11,700 

4610 

19,020 

0.68 

12,934 

17,829 

286 

0.94 

10.0 

1514.3 

56.8 

52.5 

38.0 

12,800 

5000 

21,470 

0.65 

13,956 

19,266 

311 

0.90 

11.0 

1515.8 

59.7 

55.1 

39.1 

14,000 

5450 

23,140 

0.66 

15,272 

21,055 

332 

0.91 

12.0 

1517.3 

63.0 

58.0 

40.7 

15,000 

5880 

24,880 

0.66 

16,421 

22,658 

357 

0.91 

12.5 

1518.0 

64.7 

59.6 

41.6 

15,700 

6200 

25,990 

0.66 

17,153 

23,723 

370 

0.91 
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Table 3 

Data for the analysis of heat power consumption (V160 engine, T w£ ) 625 °C). 


P m (MPa) 

n 5 (rpm) 

T c (°C) 

Nsc (-) 

Nma (-) 

Np (-) 

Co(-) 

Co/(l - T) (-) 

<Mp m V™,n s ) (-) 

3.7 

1505.0 

34.6 

7.75E+09 

0.00182 

1.41E+07 

0.2058 

0.313 

0.4225 

4.0 

1505.4 

35.8 

8.37E+09 

0.00182 

1.52E+07 

0.2052 

0.313 

0.4175 

5.0 

1506.9 

40.0 

1.05E+10 

0.00181 

1.89E+07 

0.2029 

0.311 

0.4069 

6.0 

1508.4 

43.4 

1.25E+10 

0.00180 

2.26E+07 

0.2011 

0.311 

0.4013 

7.0 

1509.9 

47.1 

1.46E+10 

0.00179 

2.62E+07 

0.1991 

0.309 

0.4007 

8.0 

1511.4 

50.1 

1.67E+10 

0.00179 

2.98E+07 

0.1975 

0.308 

0.3961 

9.0 

1512.8 

53.5 

1.87E+10 

0.00178 

3.33E+07 

0.1957 

0.308 

0.4039 

10.0 

1514.3 

56.8 

2.08E+10 

0.00177 

3.68E+07 

0.1940 

0.307 

0.3924 

11.0 

1515.8 

59.7 

2.29E+10 

0.00176 

4.04E+07 

0.1924 

0.306 

0.3894 

12.0 

1517.3 

63.0 

2.49E+10 

0.00176 

4.38E+07 

0.1907 

0.305 

0.3838 

12.5 

1518.0 

64.7 

2.59E+10 

0.00175 

4.55E+07 

0.1899 

0.304 

0.3856 



Fig. 5. Variation of electrical efficiency based on the heat power absorbed by the engine heater, as a function of the mean pressure. 


an ideal cycle without any losses caused by leakage, thermal or 
mechanical irreversibilities. This thermodynamic concept can be 
calculated for every engine if data are available about the heat 
sources’ temperatures, the dead volumes and the drive mechanism, 
according to the following functional relationship: 


the phase angle to be adjusted from 90° to 105° and consequently 
the net swept volume V ,„ to be varied from the 226 cc of early ver¬ 
sions to approximately 195 cc. 

The mechanical power losses P mec can be calculated based on 
the following model [5]: 


Co — /(Tj fc, 2 ],..., 2 (t, /d dec , fi dE , )i dR , f.i dc , /i dcc ) ( 5 ) 

The coefficients <Z> and W of indicated power losses characterise 
irreversibilities inherent to the operation at practical engine speed, 
as for example the pumping losses [4], Those factors do not depend 
on the engine speed and can be accurately calculated when the 
maximum indicated power P indi max and the corresponding velocity 
ris.max are known by means of the following expressions: 

^ _ 2(p — 3( mjx 

Nma, max 

IT / 2 w|]i!x — (q 


Notice that NMA m3X depends on the engine operating parameters 
and is inversely proportional to the indicated power loss coeffi¬ 
cients, therefore it is an index of the development level of the gas 
circuit. 

To calculate the Co values corresponding to the analysed operat¬ 
ing points, geometrical data of the V160-DMA engine are needed. 
The geometrical parameters used in the present analysis are those 
available in the PROSA 2.4 software [23], which are consistent with 
the prototype version with a modified phase angle between 
pistons whose tests were reviewed by Lista [16], This version has 
pistons with independent crank drive mechanisms, which allows 


( 6 ) 

(7) 


C mec = ae-W N * + Z + A (8) 

>3SG 

In Eqs. (3), (4), and (8), both N MA and N sc may be interpreted as 
dimensionless forms of the engine’s rotational frequency. 


4.2. Analysis of the heat power consumption 


At very low engine speeds, the indicated efficiency should tend 
to its well-known quasi-static value, namely: 


limPmd = 1 - T 


(9) 


Imposing this condition on r\ ind and taking into consideration Eq. (4) 
leads to the following restriction for the dimensionless heat power 
supplied to the engine heater: 


lim—^— 
n s ^0 p m V sw n s 


Co 

1 - T 


( 10 ) 


This constraint, not previously observed, suggests that the dimen¬ 
sionless heat power consumption QE/(p m V sw n s ) can be analysed 
based on dimensionless groups that are characteristic of the Stirling 
engine’s performance. Table 3 shows the numerical values involved 
in the analysis, which has been performed under the assumptions 
that T wE k= 625 °C and T wC « T c . Among several approaches consid¬ 
ered to match the conditions of Eq. (10), the best adjustment to 
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Table 4 

Indicated power analysis for the VI60 engine operating with helium at T wE « 625 °C. 


P m (MPa) 

n s (rpm) 

Eq. (3) (W) 

MARWEISS (W) 

PROSA (W) 

SNAP (W) 

Garcia-Granados (W) 

Error (%) 

3.7 

1505.0 

3017 

2172 

1592 

3010 

3085 

2.2 

4.0 

1505.4 

3253 






5.0 

1506.9 

4019 






6.0 

1508.4 

4780 






7.0 

1509.9 

5521 






8.0 

1511.4 

6259 

5690 

4208 

6187 

6410 

2.4 

9.0 

1512.8 

6976 






10.0 

1514.3 

7683 

6998 

5247 

7579 

7747 

0.8 

11.0 

1515.8 

8380 






12.0 

1517.3 

9060 

8403 

6134 

8896 

8880 

-2.0 

12.5 

1518.0 

9398 







the experimental data has been obtained by means of the following 
correlation, with RMSE = 0.0038 and CV(RMSE) = 4.15%: 


Qe 

p m V sw n s 


Co 

1 - T 


+ 8.871 N° u \ m N- 


( 11 ) 


Fig. 6 shows a comparison between lines of constant N MA based on 
Eq. (11) and experimental values of Q.E/(P m Vswn s ) - (o/( 1 - 1 ), 
which represents the increase in heat power consumption attrib¬ 
uted to irreversibilities. It is observed that the fluctuations of this 
value range from 7.5% to 11%, mainly due to pressure variation. 
The linear shape at practical operating conditions is caused by the 
small variation in N MA . 


are listed in Table 5. In these calculations, <P and V are assumed to 
not change meaningfully with variations of the temperature ratio t 
and the dimensionless pressure number N p , in accordance with 
previous experimental observations on very different prototypes 
[4,5], Fig. 7 also shows that the calculations based on Eq. (4) are 
in agreement with the experimental measurements of indicated 
power conducted on a recently restored V160F Stirling engine, in¬ 
stalled since the 1990s at CEDER-Soria, Spain [14]. SNAP values 
were observed to be close to the predictions of Garcia-Granados’ 
program, while those of MarWeiss and PROSA show greater 
differences. 

4.4. Analysis of mechanical power losses and alternator efficiency 


4.3. Analysis of indicated power 

Because no data were available to apply to Eqs. (6) and (7), <P 
and must be estimated from simulation software and experi¬ 
mental data. 

The results of the commercial programs MarWeiss [24], PROSA 
[23 and SNAP [25], developed decades ago, were used, along with 
results provided by the author [26) of a simulation program re¬ 
cently developed under the Envirodish Project for the SOLO VI61 
engine operating with hydrogen [15], 

Table 4 and Fig. 7 show comparisons between the different sim¬ 
ulations for mean pressures of 3.7, 8, 10 and 12 MPa, assuming 
T wE ex 625 °C for all cases. The values of indicated power calculated 
by means of Eq. (4) correspond to <P = 20 and l F = 675, which are 
the values that agree the most with the simulations provided by 
Garcia-Granados, with relatively low percentage differences, which 


Because the product of the mechanical efficiency and the alter¬ 
nator efficiency is determined after estimating the indicated 
power, obtaining each of these efficiencies can be attempted using 
Eq. (8) and data about t] a it- 

Tests of the General Motors GPU3 engine suggested expressing 
the coefficients a, [S, x, <5 of Eq. (8) by means of the following poten¬ 
tial functions [5]: 


a, x, <5 = a, 


NT, 

K 


P = 


gi 

(p l /pT 


where p and p L are the viscosities of the working gas and lubricat¬ 
ing fluid, respectively, at the reference temperature T wC i N m = mi 
RT m cl(.p m Vsw) is the dimensionless number characteristic of moving 
masses in the drive mechanism, and mi is one of those masses, used 
as a reference. 
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Fig. 6. Increase in heat power consumption due to irreversibilities. 






































































































































































































































58 


D. Garda et al. /Applied Energy 121 (2014) 51-63 



♦ Garcia-Granados (2009) 
A SNAP 
X MARWEISS 
O PROSA 

A Mesonero et al. (2012) 
-Eq.(4) predictions 


Fig. 7. Comparison of indicated power values for the VI60 engine operating with helium at T wE sa 625 °C. 


Table 5 

Mechanical efficiency analysis for the VI60 engine operating with helium at T wE sa 625 °C. 


Pm (MPa) 

n s (rpm) 

Tc (°C) 

M-) 

Nso (-) 

1 /AM-) 

Cmec ( — ) 

Pmec (W) 

P B (W) 

4 mec ( ) 

4alt ( - ) 

3.7 

1505.0 

34.6 

1.41 E+07 

7.75E+09 

1.29E-10 

0.066 

1191 

1826 

0.605 

0.684 

4.0 

1505.4 

35.8 

1.52E+07 

8.37E+09 

1.19E-10 

0.062 

1218 

2035 

0.626 

0.644 

5.0 

1506.9 

40.0 

1.89E+07 

1.05E+10 

9.56E-11 

0.054 

1308 

2711 

0.675 

0.712 

6.0 

1508.4 

43.4 

2.26E+07 

1.25E+10 

7.98E-11 

0.048 

1399 

3381 

0.707 

0.739 

7.0 

1509.9 

47.1 

2.62E+07 

1.46E+10 

6.85E-11 

0.044 

1495 

4026 

0.729 

0.785 

8.0 

1511.4 

50.1 

2.98E+07 

1.67E+10 

CD 

o 

o 

m 

1 

0.041 

1597 

4662 

0.745 

0.804 

9.0 

1512.8 

53.5 

3.33E+07 

1.87E+10 

5.33E-11 

0.039 

1706 

5270 

0.755 

0.875 

10.0 

1514.3 

56.8 

3.68E+07 

2.08E+10 

4.81E-11 

0.037 

1823 

5860 

0.763 

0.853 

11.0 

1515.8 

59.7 

4.04E+07 

2.29E+10 

4.37E-11 

0.036 

1948 

6432 

0.768 

0.847 

12.0 

1517.3 

63.0 

4.38E+07 

2.49E+10 

4.01E-11 

0.035 

2080 

6979 

0.770 

0.842 

12.5 

1518.0 

64.7 

4.55E+07 

2.59E+10 

3.85E-11 

0.035 

2149 

7248 

0.771 

0.855 


Although the engine frequency is almost a constant in these ser¬ 
ies of measurements, the mean pressure variation is sufficient to 
express the experimental data as a function of the Stirling number 
N S c = p m /(/J.n s ). This approach leads to a function for C mec deter¬ 
mined by a set of 11 coefficients, a„... , a s . Each series contains 
11 sets of measurements, but errors hinder the calculation of these 
coefficients as the solution of a system of equations. Therefore, it 
was necessary to estimate the coefficients of Eq. (8) by means of 
numerical calculus, taking into account the following information: 

- An alternator efficiency of 92% at an electrical power of 9 kWe 
has been reported for the V160F unit installed at CEDER-Soria, 
Spain, as well as efficiencies in the range of 75-80% at “low 
values” of electrical power [27], 

- The recent restoration of this unit has allowed tj a i t measure¬ 
ments to be obtained, which resolves the inaccuracy of the 
above information at low power levels [14]. 

- An efficiency of 92.5% at an electrical power of 10,850 W has 
been reported for a similar alternator coupled to the SOLO- 
VI61 unit erected under the EnviroDish Project at the CNRS- 
PROMES laboratory in Odeillo, France [20], 

Table 5 and Fig. 8 summarise the results of the best correlation 
obtained for the set of f/ oit data, which can be expressed by means 
of the following equation, assuming / i L = 0.3 Pa s: 


C mec = 0.500N p ° 13 V 7250 w w Whl)° 020 /Nsg + o.580JV p ° 199 

+ 2.99- 10 8 N p 008 /N sc (12) 

It must be noted that N m has not been considered in the list of 
variables influencing this correlation, which implies that some 
numerical coefficients are dependent not only on engine parame¬ 
ters but also on N m . However, it is thought that this simplified pro¬ 
cedure is acceptable for the objectives of this article related to the 
analysis of mechanical power losses. 

From Fig. 8 it is deduced that the correlation predicts an alter¬ 
nator efficiency of approximately 90% at an electrical power of 
7550 W, which is comparable the value of 85% estimated for a 
similar alternator on the EuroDish SOLO VI61 unit installed at 
the University of Seville, Spain, under the EnviroDish Project [15], 


5. Discussion 

Fig. 9 shows the different efficiencies of the VI60 unit as a func¬ 
tion of the characteristic dimensionless pressure number. When 
calculating the characteristic efficiency and power curves of Stir¬ 
ling engines, it is desirable to specify not only the operating vari¬ 
ables but also the parameters that remain constant, because 
hypersurfaces would be needed for a complete representation of 
models such as Eq. (3). As for the figures analysed in the previous 
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Fig. 8. Variation of the alternator efficiency as a function of the electrical power. 
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Fig. 9. VI60 efficiencies as a function of p N , operating with helium at 0.342 L x < 0.376 and N m « 0.0018. 


sections, in Fig. 9 the rotational frequency is approximately con¬ 
stant, but the temperature ratio varies almost 10% between the 
lowest- and highest-pressure points of each line. Thus those figures 
consist of projections of operating points located in different 
hyperplains. 

In addition to previous comments concerning the combustion 
subsystem efficiency and the alternator efficiency, it is observed 
that the indicated efficiency is approximately equal to 40% regard¬ 
less of the value of N p , which is influenced by the two control vari¬ 
ables (p m and T c ). It is also evident that the brake efficiency has 
small variations under medium and high pressures, where it 
reaches values slightly higher than 30%. The VI60 engine reaches 
relatively high i/ mec values compared to other prototypes, of 
approximately 77% during nominal operation in the experimental 
series examined, which has not been observed until now. 

Therefore, the mean pressure variations do not result in large 
variations in indicated efficiency, as is the case with the dimen¬ 
sionless indicated power. In contrast, the mean pressure is influen¬ 
tial in the mechanical efficiency. In agreement with Senft [28], we 
can say that the mechanical efficiency is affected by the operating 


conditions of the gas circuit, showing that forced work is greater at 
low pressures for the VI60 engine. Then the brake efficiency vari¬ 
ations with the mean pressure are caused mainly by the mechan¬ 
ical performance characteristics. 

Fig. 10 shows that efficiencies have higher variation with re¬ 
spect to the rotational frequency when the temperature ratio and 
average pressure remain constant. Despite there is one only exper¬ 
imental value of each efficiency for the considered values of r and 
N p , the curves can be obtained from the models’ equations, pro¬ 
vided that r\ ind = Cind/(Q.E/p m V 5W n s ) and rj mec = (find - £m«-)/£md. No 
experimental data are available to corroborate the curves at very 
low speeds, but the figure indicates, for example, that the engine 
could attain a brake efficiency increase of approximately 13% if it 
operates at 1000 rpm, which may be of interest from for possible 
control systems with speed regulation. 

Figs. 11 and 12 show the dimensionless power curves under the 
same conditions assumed for the above efficiency characteristic 
curves. Because dimensionless power is proportional to torque, 
this type of curves is interesting not only for power analysis but 
also from the viewpoint of operational stability. Experimental 
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Fig. 10. VI60 efficiencies as a function of Nma, operating with helium at t ps 0.374.0 and N p = 4.38 10 7 . 



Fig. 11. VI60 dimensionless powers as a function of N p , operating with helium at 0.342 ^ t < 0.376 and N MA « 0.0018. 


torque measurements of other recently developed engines [29,30] 
show trends that are consistent with the models applied in this 
article. The model of Eq. ( 8 ) is capable of explaining brake effi¬ 
ciency and torque decrease providing the mean pressure increases 
and the engine speed remains constant. Moreover, an increase in 
dimensionless indicated power with engine speed for tempera¬ 
tures and mean pressure unchanged may be observed in experi¬ 
mental engines likely due to leakage losses that occur at certain 
operating conditions. 

In summary, the models of dimensionless power and efficiency 
can be adjusted to the V160 engine experimental data with an 
acceptable degree of accuracy. These models allow us to analyse 
the influence of each operating parameter that has been selected 
as a variable, especially if the rest of the parameters remain con¬ 
stant. This procedure allows the possible margin of improvement 
for each subsystem to be evaluated. The analytical equations of 
such models are often not linear, but they have the advantages of 


being based on physical concepts and have fewer numerical coef¬ 
ficients than matrix models recently developed under the sponsor¬ 
ship of the International Energy Agency [31], 

Another advantage of using models based on dimensionless 
numbers is that they allow the results to be generalised to other 
engines using dynamic similarity criteria. 

In this way, the values of 77 * shown in Fig. 5 are clearly in agree¬ 
ment with the value of 26.5% estimated for the VI61 unit installed 
at the University of Seville, operating with hydrogen at a mean 
pressure of 11.22 MPa [15]. This result might be considered a con¬ 
sequence of the similarity between both prototypes; however, the 
VI61 engine has geometric differences, mainly at the heater, due to 
its modification for a system supplied by solar energy. Moreover, 
the use of different gases also causes the dimensionless numbers 
characteristic of the operating point to be quite different. Those 
differences reduce the dynamic similarity and are influential on 
both the indicated power and efficiency. 
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Fig. 12. VI60 dimensionless powers as a function of N MA , operating with helium at r w 374.0 and N p «4.38 ■ 10 7 . 


Table 6 

Geometric characteristics of the Stirling engines. 



V w (cm 3 ) 

vf (-) 

?*(-) 

fee-) 

Zpdx (-) 

fc(-> 


%(-) 

£(-) 

Ave (_\ 

V 2 s£ ^ 

A w r (_\ 

V 2 si 3 ' ' 

Avc 

vg* ^ 

VI60 

194.55 

0.209 

0.328 

0.190 

1.729 

0.00313 

0.00064 

0.61 

0.00313 

16.169 

972.766 

35.322 

V161 

194.55 

0.326 

0.353 

0.190 

2.005 

0.00141 

0.00093 

0.69 

0.00313 

42.038 

735.239 

35.322 


Table 6 compares the geometric characteristics of the heat 
exchangers of both engines. As stated before, the VI60 engine data 
are contained in the PROSA 2.4 software, while the VI61 data come 
from measurements made on the units installed at the University 
of Seville [26] and Odeillo [32], The 24 tubes with lengths of 
240 mm in the VI60 engine heater were replaced in the solar recei¬ 
ver of the V161 engine by 78 tubes with lengths of 320 mm and a 
smaller diameter, with the advantages of creating a heater with 
more than double the heat exchange surface and favouring a uni¬ 
form distribution of temperature, but with the disadvantages of 
increasing the dead volume and decreasing the hydraulic radius- 
to-length ratio. 

An indicated power of 13,670 W and indicated efficiency of 48% 
were calculated for the VI61 unit of the EnviroDish Project, oper¬ 
ating with hydrogen at p m = 11.22 MPa and T wE ~ 732 °C [15], It 
must be noted that the authors question the validity of this 
mechanical efficiency, which should be equal to 64.9% to attain 
the measured electric power of 7550 W with the estimated alter¬ 
nator efficiency of 85%. 

As discussed below, this mechanical performance actually 
seems low because the VI60 and VI61 engines have similar drive 
mechanisms and Eq. (12) predicts that the V161 engine mechanical 
power losses would be of the same order as the ones for the V160 
at high pressures (Fig. 13). Furthermore, the inconsistency is even 
greater when the higher alternator efficiency for the electric power 
measured predicted by Fig. 8, of approximately 90%, is considered. 

It is also interesting to note that the lines in Fig. 13 have the 
appearance of Stribeck’s curves, where friction increase at high 
pressures and low velocities may be caused by changes in the 
lubrication regime. It is also evident the analogy between the Her- 
sey number and the inverse Stirling number, which has already 


been argued as a foundation of the model of Eq. (8) [5,9], Although 
there are no experimental evidences for values greater than 
12.5 MPa, the trend lines of Figs. 5 and 9 indicate the possibility 
that the mechanical efficiency will decrease above a certain mean 
pressure value, consequently affecting the characteristic curves of 
brake and electrical efficiency. Similarly, the model of Eq. (8) 
may explain the decrease in brake power at relatively high mean 
pressures observed in other recent experiments already cited 30], 

Eq. (4) provides a way to analyse whether the VI61 engine can 
reach the indicated power of 13,670 W predicted by the simulation 
program, and then assess whether the mechanical efficiency 
should be higher. 

First, with the geometrical data in Table 6 and the parameters 
corresponding to the operating point considered, listed in Table 7, 
the value of Co that is obtained for the VI61 engine is greater than 
those obtained for the V160 at high pressure levels. 

Furthermore, it is expected that the loss coefficients <P and W 
are lower for the VI61 engine operating with hydrogen than for 
the VI60 engine with helium because the following relationship 
was observed for the GPU-3 engine operating with both gases [7]: 

jlli H 2) -.46 q 92 

0 . f'(He) ~ 50 

In addition, the differences between the and coefficients of 
the VI60 and VI61 engines can be justified by geometrical reasons. 
As shown in Table 6, the ratio r hR /L R and the regenerator porosity 
are 45.3% and 13.1% higher, respectively, in the V161 engine than 
the assumed values for the VI60. Because the indicated power 
losses usually are an order of magnitude higher in the regenerator 
than in the other heat exchangers [33], it is expected that the indi¬ 
cated power loss coefficients of the V161 engine are lower than 
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INVERSE STIRLING NUMBER 


Fig. 13. Variation of dimensionless mechanical power losses as a function of the inverse N S g with N p as a parameter. 


Table 7 

Indicated power analysis for the VI61 engine operating with hydrogen at 
T w e ~ 732 °C. 


Pm (MPa) 

n s (rpm) 

Nma 

Nsg 

N P 

T c (°C) 

Co 

11.22 

1523 

0.00128 

4.61 E+10 

5.87E+07 

55.3 

0.2019 


those of the VI60 engine, even operating with helium. However, 
even with zero loss coefficients, 13,670 W cannot be reached using 
Eq. (4). 

Whatever the origin of the discrepancies, it can be seen that Eq. 
(10) leads to P mec rs 1970 W for the above operating point of the 
V161 engine, as well as P B ~ 8389 W, assuming 90%. Thus, 
Pin d rs 10,359 W and t7 mec ~81%. 

Eq. (4) is in agreement with this indicated power calculation 
using <P rs 13 and 'F rr 440, which leads to rs 0.0060. This 

value of N M A,max exceeds the results observed for other engines un¬ 
til now, which indicates than the VI61 gas circuit is at an advanced 
developmental level, even if the indicated power is well below the 
predictions of the simulation program. Higher loss coefficients 
would result in a smaller indicated power, but in this case the 
mechanical efficiency or the alternator efficiency would have to 
be even greater. 

The possibility of achieving Nma, max ~ 0.0060 can be evaluated 
by experiments carried out with the EuroDish unit in Odeillo be¬ 
cause the <P and '¥ coefficients vary little with changes in temper¬ 
ature and mean pressure. In this case, P e a;10,850W was 
measured at a moderate direct normal irradiance level of 
906 W/m 2 and low ambient temperature of -5 °C [20]. The engine 
operating point is not reported by the authors, but the values of 
T wE rs 780 °C,Twc ~ 40 °C and p m rs 14 MPa are deduced from an¬ 
other publication related to the same experiments [34], Substitut¬ 
ing the values N S c ~ 5-72 ■ 10 10 and N p rs 7.29 ■ 10 7 into Eq. (12) 
allows P mec r«2517W to be obtained for these operating condi¬ 
tions. Therefore, the assumption of ri a it~ 92.5% leads to 
Pm ~ 14.247 W, and consequently ii mec rs 82%. 

A quasi-static dimensionless indicated power of Co ~ 0.2172 can 
also be computed for the operating conditions of the Odeillo unit, 
so that Eq. (4) allows P, nd rs 13,884 W to be obtained at n s - 
rs 1530 rpm for the coefficients of indicated power loss estimated 


at the Seville unit, i.e., cP«13 and W rs 440. The difference of 
2 .6% between the indicated power calculations may be due to 
small errors in the measurements of operating conditions or, more 
most likely, in the estimations of <P or '¥. However, it can be 
concluded that the estimates should be close to reality and the 
Odeillo data are consistent with a value of N M A,max slightly greater 
than 0.0060. 

Therefore, the indicated power of the VI61 engine may exceed 
the value of 13,670 W predicted by Garcia-Granados’ simulation 
program, but not at the operating conditions considered for the 
unit installed in Seville. It seems unlikely that the discrepancy is 
due to a fundamental flaw of the simulation program, as in previ¬ 
ous sections it has been seen that its predictions proved acceptable 
for the VI60 engine, but the simulations may have been performed 
with inadequate assumptions about the geometric characteristics. 
For example, the authors mention [15] that they adopted geomet¬ 
rical data specified by Organ 35] for the V160 engine. Those data 
refer to an engine version with a heater dead volume somewhat 
lower than the value specified in Table 6 for the VI60 engine, 
and, more importantly, with a drive mechanism that has a phase 
angle of 90°, and therefore V* rs 225 cm 3 and different instanta¬ 
neous volumes. 

In short, the analyses show that both prototypes have an 
advanced degree of development as evidenced by relatively high 
efficiencies for each subsystem. 

6. Conclusions 

VI60. engine tests have been performed with small variations 
in rotational frequency, but with mean pressure variations suffi¬ 
cient to achieve the research objective of deducing power and per¬ 
formance characteristics for the internal subsystems. 

The combustion subsystem shows very uniform efficiencies for 
all series of measurements, with values of approximately 90%. 

The internal subsystems have been analysed by means of simi¬ 
larity models which have been previously developed through 
experimental data from very different prototypes and are consis¬ 
tent with recent measurements in other experimental engines. 
Thermodynamic models for the gas circuit are also taken into ac¬ 
count. Model uncertainties have been evaluated by considering 
experimental data for indicated power and alternator efficiency. 
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For the VI60 engine operating with helium at a mean pressure 
range from 10 to 12.5 MPa, values of electrical efficiency, based on 
the heat power absorbed by the engine heater, are slightly higher 
than 26%, which is in agreement with calculations for the VI61 
unit installed at the University of Seville, operating with hydrogen 
at a mean pressure of 11.22 MPa. 

The increase in heat power consumption which is attributed to 
irreversibilities, varies from 7.5% to 11%, mainly due to the varia¬ 
tion in pressure. The curves of heat power absorbed by the engine 
heater as a function of the mean pressure show a linear shape at 
practical operating points due to the small variation of the charac¬ 
teristic Mach number. 

Operating with helium at a mean pressure of 12 MPa and a tem¬ 
perature ratio of 0.374, the VI60 engine can achieve an indicated 
efficiency of 40%, with a mechanical efficiency of 77%. Therefore, 
it produces 5875 W of electrical power at approximately 
1517 rpm, with an alternator efficiency of approximately 84%. 

The electric power of 7550 W measured at 1523 rpm on the 
VI61 engine installed in Seville, operating with hydrogen at a 
mean pressure of 11.22 MPa and temperature ratio of 0.327, is 
compatible with a mechanical efficiency of 81% and an alternator 
efficiency of 90%. These values correspond to an indicated power 
of approximately 10,360 W, significantly lower than that estimated 
in previous publications, most likely due to errors in the geometric 
parameters assumed. 

The electric power of 10,850 W measured at 1530 rpm on the 
VI61 engine installed in Odeillo, operating with hydrogen at a 
mean pressure of 14 MPa and temperature ratio of 0.297, is 
compatible with a mechanical efficiency of 82% and an alternator 
efficiency of 92.5%. 

Therefore, the analyses show that both prototypes have an 
advanced degree of development as evidenced by the similar, 
relatively high efficiencies of each subsystem. 

Nevertheless, the engines are not dynamically similar. From this 
point of view, the differences between the VI60 and VI61 engines 
are not important in terms of geometric parameters, temperatures 
or pressures, but it is remarkable how much is influenced by the 
type of gas. Consequently, the indicated power losses caused by 
irreversibilities are greater in the VI60 engine than in the VI61, 
whose operating point corresponds to a characteristic Mach num¬ 
ber a third smaller. In addition, the VI61 engine coefficients of 
indicated power losses are so low that the characteristic Mach 
number corresponding to the maximum indicated power point 
reaches one of the highest values observed so far (N MAmax « 
0.0060). On the other hand, the use of hydrogen increases the 
characteristic Stirling number of the VI61 engine, which allows 
moderate mechanical losses to be achieved at pressures even 
higher than those of the V160 engine. 
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